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SUMMARY 


The potential curves of the lowest states of HZ have been calculated in the range of R-values 
from 1 to 10 a.u. The following states, arranged according to increasing value of the minimum 


_ energy, are found to correspond to discrete energy levels: *D}, *I1,,*X,#X* and4X}. The minimum 


oJ 


of the potential curve for the ground state is found at R=1.65 a.u. and E = —1.1057 a.u., corre- 
sponding to a dissociation energy D,=59 kcal/mole. Hz is thus found to be comparatively stable 
towards dissociation into H and H-. Moreover, the electron affinity of H, is found to be EA(H,) = 
—0.28 ev = —6.4 kcal/mole. The bearing of these results on the possible reactions between elec- 


‘ trons and H, molecules is briefly discussed. 


1. Introduction 


In a recent discussion of certain dark~structures in the solar corona Ohman [1] 
has suggested that these structures may be explained as absorption structures, 
originating from H- ions. According to Ohman one possibility for the formation of 
these ions would be by collisions between electrons and H, molecules in the ground 
state, followed by capture of the electrons in formation of Hz ions with subsequent 
dissociation of Hz into H and H-. 

Starting from a suggestion by Ohman that inverse predissociation should be 
important for the recombination of H,, the present author [2] has suggested some 
alternative processes for the formation of H~ by way of Hz. For a more detailed 
study of the transition probabilities in these reactions the potential curves of the 
ground state and the lowest excited states of Hz must be known with some accuracy. 
The present investigation aims at a theoretical determination of these potential 
curves. 

For a long time the existence of the negative hydrogen molecule ion has been 
rather questionable, as Hz had not been observed in any experiment. However, a 
recent investigation by Khvostenko and Dukelsky [3] of formation of H- ions from 
collisions between electrons and hydrogen molecules indicates that Hz ions are possibly 
formed as an intermediate step in the reaction. Moreover, in a very recent extension 
of this work [4] the authors claim to have ascertained the appearance in the mass 
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spectrometer of Hz in a yield 20-30 times lower than the yield of H-. These experi- 
ments will be discussed further in connexion with the results of the present investiga- 
tion. 

Previous theoretical investigations of the Hz ion by valence bond (VB) methods 
have been carried out by Eyring, Hirschfelder and Taylor [5] and recently by Dal- 
garno and McDowell [6]. Eyring et al. [5] have studied the lowest state, 231, OF Ua 
in connexion with an investigation of the ortho-para hydrogen conversion by alpha- 
particles. By use of a VB treatment with two variation parameters they found that 
the ground state of Hz should be stable, although the formation of the ion should 
have low probability. They obtained a potential curve with a minimum at an inter- 
nuclear distance of 1.8 A, the minimum being 20.9 kcal/mole lower than the potential 
computed for infinite separation of the nuclei. This limiting value, corresponding to 
H +H, was however 37.4 kcal/mole higher than the energy of H + H + e, implying 
that the H- ion should be unstable in contradiction to the well-known fact that the 
ground state of H- is stable, having a dissociation energy of 17.40 kcal/mole [7]. 
Therefore, the potential curve of Hz, computed by Eyring et al. [5], is not so well 
suited for a closer analysis of reactions involving Hz ions as well as H, molecules 
and free electrons. For this reason the present investigation is concerned not only 
with excited states of Hz but also with a reexamination of the ground state. 

The recent investigation by Dalgarno and McDowell [6], concerned with the 
ground state and the lowest 2X5 state of Hz, is more elaborate than [5]. Still it takes 
only a single configuration into account. It was published, while the present study 
was in progress, but unfortunately not noticed by the author until her own calcula- 
tions were almost completed. Although the present treatment has been carried out 
independently of [6] the first part of it follows essentially the same lines as [6] and 
gives almost equivalent results. These results will nevertheless be presented here, 
because they will be given in a more informative manner and because they are 
essential for the following part, which takes into account superposition of configura- 
tions. A detailed comparison of the methods used in [6] and used here will be given 
below. 

In Section 2 of the present paper an examination of the possible lowest states of 
Hz is given. In Sections 3 and 4 appropriate electronic wave functions for the ground 
states of H- and Hz are discussed. In Section 5 the wave functions for some of the 
lowest states of Hz are given. Section 6 presents the results of the calculation of the 
potential curve for the ground state, when a single configuration is used. Section 7 
contains the results for the 2X, states obtained from superposition of configurations. 
Section 8 contains the corresponding results for the 2X; states. In Section 9 the results 
for the *II,,, states and for the quartets are given. A discussion of the results and of 
the bearing on the possible reactions between electrons and hydrogen molecules is 
given in Section 10. The Appendix contains some remarks about the calculations 
of the pertinent integrals and gives the computed energy values. 


2. Classification of the lowest states of Hy 


In order to find the potential curves of the lowest energy states of the Hz ion 
it is suitable to start from the lowest states of the dissociation products H and H-. 
In fact, the large equilibrium distance of Hy, found by Eyring et al. [5], and the still 
large distance, found by Dalgarno et al. [6], indicates that the states of the molecu- 
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lar ion should more likely be similar to the states of the separated atoms than to the 


states of the united atom. 

Now, as to H-, only the ground state has been investigated previously. This 
state is a 1s? 148 state. The lowest excited states of heliumlike atoms are 1s2s 319 
and 1s2p*"P. As is well known the energies of the different states of two-electron 


_ atoms are approximately proportional to (Z — s)?, where Z is the nuclear charge and 


: 
4 
, 
a 
i 
if — 
oF 


4 
if 
= 
4 


s the screening constant. Thus, it can be expected that the lowest excited states of 
H- should be very close to the ground state. However, as H~ has a rather small 
electron detachment energy, it could well be that all excited states of H- belong to the 
continuum. Nevertheless, it has been assumed that the lowest states of Hz should 
be derived from the ground state 1s 2S of H combined with the ground state and the 
lowest excited states of H-. : 

The probable sequence of the lowest states of H~ can be obtained in the following 


_ way. As is well known from experiments, the energy of the ground state of iso- 


electronic atoms is to a remarkably good approximation given by the formula 
H=AZ+B2Z+C. (1) 


In fact, when A, B, C are determined by means of the ground states of He, Lit, 
Be*?2, B*3, C+4, it is e.g. found that the ground state energy of H~ should be 16.91 
kcal/mole, in excellent agreement with the value 17.40 kcal/mole of Hylleraas et al. 
[7]. Assuming the formula (1) to be valid also for corresponding excited states and 
by application to the above-mentioned series of atoms it is found that the sequence 


_ of the excited H= states presumably is the one, given in Table 1. It should be men- 
- tioned, however, that the accuracy of the values is not good enough to ascertain 


the actual existence of any discrete excited states of H-. 

The states of the molecular ion Hz, which by group-theoretical methods are ob- 
tainable from the given atomic states, are also collected in Table 1. It should be 
remembered, however, that all these states do not necessarily correspond to actual, 
descrete molecular states, as some of them may belong to the continuum. 


3. Ground state of H™ 


The ground state of H~ has been thoroughly investigated by several authors. The 
object of the present. study is not to achieve a more accurate value of this ground 
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~ 7 a = = 
state energy. What we are looking for is a wave function, simple enough to be used 
as a basis for the wave functions of Hz but accurate enough to show the stability 
of H-. : 
The simplest wave function for the 1s? 1S state is 


@, = Det [xa (¢; 1)%a(2;2)|. . (2) 


— 


Here, the notation Det|...| stands for the determinant, having the diagonal term .... 
The bar denotes an orbital with negative spin. Absence of the bar implies positive 
spin. 7.(¢;i) is an 1s atomic orbital (AO) of the electron 7 at the atom a and with 
the orbital exponent C: 


=) ~ 
(tsi) |E exp (—E ra) (3) 


In the following, the arguments of the AO’s will not be written down explicitly: 
but in any product the arguments of the AO’s will always be in the sequence 1,2,3.... 

The function (2) was used by Eyring et al. [5] for the limit of infinite separation 
between the nuclei of Hz. As mentioned above, they did not obtain stability for H—- 
by this function, although they treated ¢ as a variation parameter. Now, it is well 
known from the study of heliumlike atoms by e.g. Hylleraas [7] that the computed 
energy value can be considerably improved by choosing different orbital exponents 
for the two 1s AOQ’s: 


®, a Det | 7a (C) Za(C")| 7 Det| a (C’) Ha(C) | : (4) 


This function has previously been used for H- by Chandrasekhar [8]. 

Treating ¢ and ¢’ as independent variation parameters, he found a value of the 
energy of H-, lower than the energy of H + e, i.e. H~ is found to be stable. Moreover, 
he found ¢’ to be much smaller than ¢, about }¢. This indicates that the single con- 
figuration 1s? does not describe the 1S state of H~ very well. Consequently, if we had 
chosen to write the wave function of the 4S ground state as a superposition of con- 
figurations, composed of hydrogen AO’s, we would have found a large contribution 
from the configuration 1s2s. This conclusion is in agreement with the calculation 
on He by Taylor and Parr [9]. By inspection of Table 1 it is seen that the 4S ground 
state of H~, as described by the wave function (4), can be interpreted as a superposi- 
tion of the two configurations in Table 1, which give rise of 18 states. Obviously, 
these two configurations will also give rise to an excited 18 state, orthogonal to the 
ground state. It can however be predicted that this excited 18 must have a rather 
high energy value and thus will not contribute essentially to the lowest ¥ states 
of Hz. Consequently, the excited 4S state will not be considered in the following. 


4. Ground state of Hy 


The Hz ion is a system composed of two nuclei a and 6, separated by a distance R, 
and three electrons 1, 2 and 3. We are looking fora comparatively simple eigen-func- 
tion of the hamiltonian H of this system. As we are interested in the whole potential 
curve, the VB method should be the most suitable one. According to this method in 
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e event form the unnormalized, electronic wave function of the 2Xi ground 
state is: 


¥F, = Det | 7a(2)Za(0) x0 (6)| — Det | yo (0) Z0()xa(C)|, (5) 


where (C) is the function (3). It is easily seen that in the limit R =o ‘Y’, takes the 
‘form of ®,, multiplied by the 1s function of the free H atom. 
It is of some interest to compare the VB function (5) with the function, obtained 


a the simplest molecular orbital (MO) method. From 1s AO’s two MO’s can be 
_ formed: 


= hoe 


;. Yr (C) =Xa(C) + Xo (6), 6) 
. | yo (£) =a (E)— to (C). 

_ Here y, is a bonding MO and y, an anti-bonding MO. If we choose equal orbital 
i exponents of y, and ws, the total wave function of the ground state can be written: 


Fy = Det |p, P,y2| = Det] (ya + 40) (Za + Zo) (Ha — Xo) | = — 24. 


_ Thus, it is seen that, disregarding a numerical factor, which will disappear by proper 
normalization, the VB method and the MO method in their simplest forms will 

_ give the same ground state wave function of Hz. This function is also identical with 
the wave function, used by Eyring et al. [5], and with the simplest of the three 
different functions, used by Dalgarno and McDowell [6]. 


JWT EV, de 
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The electronic energy 


(7) 


is a function of the internuclear distance # and of the orbital exponent ¢. The simplest 
way of calculating the potential curve is to choose the ¢-value for the free atom and 
calculate H# for different R-values. A refinement can be obtained by variation of ¢ 
for each chosen R-value. This refined calculation has been carried out previously 
by Eyring et al. [5] and by Dalgarno and McDowell [6] but has been repeated by the 
present author to make possible a detailed comparison of the results of this treatment 
with the results of alternative calculations.. 

As has been pointed out above, the limiting function ®, of the function ‘f’, does 
not describe the stability of the H— ion correctly, not even when ¢ is varied. To 
obtain a more correct limiting value of the potential curve we must choose a wave 
function of Hs, which takes the form of ®, in the limit R = oo. This requirement is 
fulfilled by the following function: 


Ys = 1 —%2 + Ps — Pa (8) 

where p, = Det | a (5) Za (0) 40 (C") | 
G2 = Det | %p (2) Zo (6') Xa (6) |, | (9) 

g3= Det Oke ee ane | 

pa= Det |%o(6') Zo (6) Xa (5"") |. J 
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_ The electronic energy can now be expressed by means of the matrix elements of H, 


Hy={ oi H oat, (10) 
and of the overlap matrix, 
An={ of Gr adt. (11) 


The energy will be a function of the four parameters R, ¢, ¢’ and ¢’’. To obtain the 
potential curve the three orbital exponents have been varied independently at each 
chosen R-value. 

The function (8) is equal to the best wave function, used by Dalgarno and McDowell 
[6]. However, these authors did not vary all the orbital exponents independently. 
Instead they chose the exponents to be equal to the values, obtained for Rk =o, 
multiplied by a common scale factor 9, which they varied for each chosen R-value. 


5. Lowest excited states of Hy 


As is shown in Table 1, Hz states of the same symmetry can be obtained from 
several states of the separated atoms. As the states of H- are expected to lie very 
close together, it can be anticipated that Hz states of equal symmetry will interact 
considerably. 

In principle, wave functions of states with different symmetry could be written 
down independently, and the variation parameters could be chosen differently. How- 
ever, in order to keep the computational work within tolerable limits, we have 
proceeded along the following lines. 

From the computation of the potential curve of the ground state we have found 
that for all R-values the orbital exponents ¢ and £”’ are both almost equal to 1.00. 
Moreover, the remaining orbital exponent ¢’ varies slightly with R but is rather 
close to 0.30 in the whole interval investigated (R =4 a.u. — co). Extrapolating this 
result to the calculation of the excited states, we will not treat the orbital exponents 
as variation parameters, but give them the values 


C=" =1.00, ¢’ = 0.30 (12) 


both for different R-values and for states of different symmetries. Instead of the 
orbital exponents the variation parameters will be the coefficients of the different 
configurations, taken into account. 

The remaining question is how to choose the orbital exponent ¢’”’ of the 2p AO of 
H-. This 2p AO should describe the polarization of H- in the field of the H atom. 
Obviously, the very diffuse AO 1s(£’) will be polarized much more easily than the 
concentrated AO 1s(¢). Consequently, the 2p AO should belong to the set of diffuse 


AO’s, i.e. €’”’ should be 
= 0.15. (13) 


The choice (13) should be the best at least for the 2, ground state of H-. 
Using this value we find the energy of the 1s2p °P state of H- to be —0.4891 a.u. 
This value could be compared with the value obtained from extrapolation by aid 
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C5 /a 14 COS Og exp (—(0" Fe); 


——___ ————_ oe ‘ 
Ta=VC""*/a rq sin A, Cos Gq exp (—C'” 14). J 


€ ee it will be assumed that the positive lobes of o, and a, are directed towards 
€ AL otner. ate 
_ Using the notation (14), the configurations, essentially contributing to the lowest 
_ states of Hz, can be written: 


~ 


2 = Det | 5.34 s9|, py= Det | 5), 35 8s, 
P2= Det | 5,85 |, @,= Det | 84 8.30... (15) a 
C= Det | 858.80 |, Pg = Det | 85 85 Sa\; 

¥ G7 = Det | 82 G 85\, Pig = Det | oy 5p Sal, 

E Gg = Det | 8 Gy 8a; P11 = Det | Sa 6a So, (16) 

4 Pq = Det | oa 5a 8p |, P12 = Det | 5p on Sa); 

g 

p Gis— Det |sadasl; Ge = Det |r035 5a, 

Gia = Det | 8 Zp 8a; G17 = Det | 8270030, (17) 

J P15 = Det | 05a 80, Pig = Det | 5, 2 Sal. 


It is easily seen that 91, Po, 3, Ps Of (15) are identical with those obtained from 
(9) after introduction of (12). In the limit R =o these functions are obviously con- 
nected with the function (4), describing the ground state 4S of H-. 9, Pe, Ps Pa will 
therefore describe the interaction of an H atom with an H~ ion in its stable 1S 
state. y; and gy, on the other hand, will describe the interaction of an H atom with 
an H- ion in a hypothetical #8 state. Similarly 77, Ps, Po. Pio OF (16) as well as G5, 
— Pra Yrs» Pig Of (17) will describe the interaction of an H atom with an H~ ion in a 
_ 1P state. Finally 91, Py, of (16) as well as gy, Pig of (17) will describe the interaction 


of an H atom with an H~ ion in a ?P state. 
To find the correct combinations of these different functions, which actually 


describe the molecular states: of interest, we let the spin operator S* operate on the 
different y’s. In this way we obtain the following wave functions for the different 


states of Hz: 
39 
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V3 = 91+ Pst (Pat Pa); 
a 1, 
pie + (H Ye (18) 
“tes = 7+ Pot (Pst Pio)» 
— Py + (Ps — Pia) 
Pee — 3+ Ps + (Y2— Pat Pe)s 


Br (19) 
— Pot Put (Ys — P10 + Piz); 
TL, _ = Pist Pris = (Pra + Pie): (20) 
P10 = Pis — Piz £ (Pre — Fis); 
Iu. Cia = P13 — Pis + Piz + (Pia — Pie + Pis)- (21) 


Here the upper signs give the w states and the lower signs the g states. It is seen that 
(18)-(21) will give all the states of Table 1, derivable from the H- 18 ground state 
and the hypothetical H- states ?P, 1P and 3S. 

‘To obtain the best description of the lowest state of each symmetry, linear com- 
binations have been made of all functions of the same kind. As mentioned above, 
the coefficients in these combinations have been chosen as variation parameters and 
determined by minimizing the energy. 


6. Results for the ground state by use of a single configuration 


In this Section we will discuss the results of the calculations on the ground state 
according to Section 4. The calculations have bee carried out in four different approxi- 
mations: 


I the same value ¢ of all orbital exponents, this value being fixed to the value 

for the free H atom; 

II the same value ¢ of all orbital exponents, this value being determined by varia- 
tion for each internuclear distance; 

III two values ¢’ and ¢=¢"' of the orbital exponents, these values being determined 
by variation for mihi internuclear distance; 

IV three values ¢, ¢’ and ¢”’ of the orbital exponents, these values being determined 
by variation for each internuclear distance. 


I and II have been computed by use of the electronic wave function Y, of (5), and 
III and IV by use of the wave function ‘’; of (8), where in case III we have put 


pee - 


ie eae 

II is identical with the treatment of Eyring et al. [5] and with one of the cases 
considered by Dalgarno and McDowell [6]. These authors have also calculated a 
case identical with the present case I. Case IV is similar to the most accurate treat- 
ment in [6], although in [6] all three orbital exponents were assumed to be of the form 
¢ =¢.0, 80 that only the common parameter g had to be determined by variation 
for each internuclear distance. 

The results of the energy computation are shown in Fig. 1, where the calculated 
energies H are plotted against the internuclear distance R. The horizontal lines 
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Fig. 1. Potential curves for the ground state 2X; of Hy, obtained from the single configuration 


(8). I: €=C’ =C” =the value for the free H atom. II: Variation of one exponent: €=¢’=C”. III: 

Independent variation of two exponents: ¢’ and =’. IV: Independent variation of three 

exponents: ¢, ¢’ and ¢”’. V: Value of H + H~- according to [7]. D = Dissociation energy; the lengths 

of the vertical arrows show the computed values of D; the positions of the arrows indicate the 
computed equilibrium distances. 


indicate the energy values at infinite nuclear separation, H,,, in the different approxi- 
mations. The lengths of the vertical arrows, D,, Dy, etc., indicate the computed 
values of the dissociation energy. The positions of the arrows show the computed 
equilibrium distances. The pertinent numerical values are listed in Table 2 in the 
Appendix. 

Fig. 1 shows that the computed values of H#,, are improved in the sequence I, 
II, III, IV. These values should be compared with the energy level V, which is the 
sum of the best energy value of H-, computed by Hylleraas et al. [7], and the energy 
value of the H atom. It is seen that the introduction of one variable orbital exponent 
improves the value of Z,, by as much as 0.0654 a.u. = 41.1 kcal/mole. The introduc- 
tion of two variable orbital exponents improves this energy value by another amount 
of 0.0722 a.u. = 45.3 kcal/mole, giving an energy of H + H- well below the energy 
—1.00 a.u. of H+H-+e in agreement with the well-known stability of H~-. The 
introduction of three independent orbital exponents gives a further improvement 
of the energy value by only 0.0007 a.u. = 0.4 kcal/mole. 

Furthermore, Fig. 1 shows that the values of the computed dissociation energies 
are in the order D, > D,, > Dy, > Dry. The differences between D,;, Dy and Dy; 
are considerable. Thus, the potential curve, obtained by use of the simplest wave 
function, becomes remarkably modified, when the wave function is improved. More- 
over, it is seen that D,y is almost equal to D,1, indicating that no serious error should 
be introduced by putting ¢=¢” as has been done in Section 5. 
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Fig. 2. Values of orbital exponents ¢. I: Constant value of €=C’=C”’ =the value for the free atom. 
II: Variation of one exponent: €=C’ =”. III: Independent variation of two exponents: C=C” 
and ¢’, IV: Independent variation of three exponents: ¢, ¢’ and ¢”. 


It should be noted that although the absolute energy values will always be im- 
proved when the wave function is improved, this is not at all true for the computed 
D-value. As long as the value of #,, is far from the correct one, the D-value may 
well diminish drastically, when the wave function is improved. 

In Fig. 2 the values of the orbital exponents, obtained in the four different approxi- 
mations, are shown in terms of R. It is seen that the mean value of ¢€, ¢’ and €” is 
about the same for all the three cases II, III and IV. This could be helpful for finding 
the minimizing values, when more orbital exponents are varied independently. 

Curve IT of Fig. 1 is the same as that published in [5]. In [6] the authors have only 
published the energy of interaction between H and H-, i.e. H — E#,,. At least as far 
as this difference is concerned the curves I and II of Fig. 1 seem to coincide with 
the corresponding curves of [6]. Moreover, curve III of Fig. 1 is very nearly the same 
as the curve obtained from the most accurate treatment in [6], although the two 
approaches are slightly different. The large excess labour, involved in the present 
approach IV, apparently does not pay, as it hardly changes the D-value at all. 

According to the results given in this Section, the energy of interaction between 
H and H~ should be of the same magnitude as van der Waals’ energies. If these 
results were conclusive, Hy could hardly be considered as a stable molecule and its 
importance for the formation of H~, as suggested in [2], would be rather questionable. 


7. Results for the ground state and other 2X, states from superposition 
of configurations 


In this Section the results of the calculations on the 2X; states according to Section 
5 will be presented. 
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_ Fig. 3. Potential curves for the ground state *X of H,, obtained from superposition of configura- 


tions, given in (18). A —— Energy obtained from ¥’;, composed of 1s AO’s only. B —:— Energy 

_ obtained from superposition of Y, and ‘’,, composed of 1s AO’s only. C, D === Energy obtained 

from superposition of Y’;, Y,, Y; and 4, composed of 1s and 2p AO’s. H, Energy of 
H,, computed by the VB method. 


2 4 6 8 


The calculations on the ground state have been carried out in four steps by use 
of an increasing number of the wave functions, given in (18): 


(A) Y, only; 2 
(B) superposition of ’, and ’,; 

(C) superposition of 3, ‘, and V’,; 

(D) superposition of ‘¥,, ‘,, ‘, and ‘Ys. 

The first step A is almost identical with case III of Section 6, although the constant 
values (12) of the orbital exponents are used in case A. In Fig. 3 the result of this 
calculation is given by curve A of dashes. A comparison of this curve with the curve 
III of Fig. 1, or of the energy values, listed in Tables 2 and 3 in the Appendix, shows 
that the potential curve A is very nearly the same as the curve III. 

The step B consists of a superposition of the functions ‘’; and ‘",, both composed 
of 1s AO’s only. This superposition describes the interaction between the 727 
states, derivable from the H- ground state 18 and from a H~ 3S state. In Fig. 3 
the result of this superposition is given by curve B of dots and dashes. For values 
of R>5 a.u. the contribution of , to the ground state is almost negligible. At 
smaller R-values however the superposition changes the potential curve drastically, 
giving rise to a new and much lower minimum at & = 1.78 a.u. i 

The step C consists of a superposition of ‘Y’;, ‘i’, and ‘’,, describing the interaction 
between the 2X; states, derivable from the H- 1S, °S and *P states. The step D, 
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E a.u. 


Fig. 4. Potential curves for the lowest #27 state of H,, obtained from superposition of con- 


figurations, given in (18). A — — Energy obtained from ‘’;, composed of 1s AO’s only. B —--— 

Energy obtained from superposition of Y’, and ¥’,, composed of 1s AO’s only. C == Energy 

obtained from superposition of Y’,, Y", and ‘’,, composed of 1s and 2p AO’s. H, —— Energy 
of H,, computed by the VB method. 


finally, consists of a superposition of ‘Ys, ‘Y",, ‘; and Y’,, describing the interaction 
between the 2; states, derivable from the H- 1S, 3S, 1P and °P states. The results 
of the two last steps are so similar that they are both represented by the thick, 
solid curve in Fig. 3. The pertinent numerical values are collected in Table 3 of the 
Appendix. The thin, solid curve in Fig. 3, as well as in Figs. 4, 5 and 6 below, repre- 
sents the potential curve of the ground state X1X; of H, and is only given for the 
sake of comparison. For this reason the curve has been computed by a VB treatment 
of an accuracy, corresponding to the accuracy of the Hz calculations, i.e. by use 
of 1s AO’s with a constant value 1.00 of the orbital exponent. No 2p AO has been 
included in the H, calculation, as the 2p AO of Hz only has reference to the polariza- 
tion of the outer, diffuse 1s(¢’) AO. 

The inclusion of 2p AO’s in the wave function for the ground state of Hz is seen to 
deepen the minimum considerably and to shift it somewhat towards smaller R-values. 
The minimum is now found at Rk = 1.65 a.u. and # = — 1.1057 a.u. It is found to 
lie very close to the computed minimum of H, and only 0.0102 a.u. = 6.4 kcal/mole 
above this minimum. The D,-value, computed for Hz according to the final step D, 
is 0.0933 a.u. = 58.6 kcal/mole. Thus, the present calculation shows that Hs should 
be about as stable towards dissocitaion as molecules in general. On the other hand, 
the computed value of the electron affinity EA of H, is found to be negative, although 
it has been reduced by a factor ten in comparison to the results of [6]. The value 
according to step D is EA(H,) = — 0.28 ev = — 6.4 kcal/mole. 
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i end I 4 - are il * 
r differs nt wave functions of *L; states, given in (18), can, of course, 
e used for the computation of potential curves of three excited 2D} states. 
lowever, the different wave functions of (18) are not orthogonal, so that the poten- 
al curves, obtained from the wave functions V,, V';, V, separately, will have no 
physical significance. A step, corresponding to step A, can therefore not be used for 
the calculation of any excited state. On the other hand, in the steps B, C and D one 
_ obtains orthogonal, linear combinations of Y’,, ¥',, Y, and W,, two in step B, three 
_ instep C and four in step D. In these steps potential curves can therefore be calculated 
for one, two and three excited states respectively. 
. The energy values of the three excited 2X; states, obtained from step D, are listed 
_ in Table 4 of the Appendix. Two of the states are repulsive. The third excited state 
_ has a minimum at & = 1.73 a.u. and an energy value of — 1.0472 a.u., corresponding 
— to D, = 0.0581 a.u. = 36.5 kcal/mole. 
_ Moreover, it is of some interest to study the coefficients of the different ¥’s in 
__ the linear combinations and their variation with R. It has been found that for R-val- 
ues of 1-2 a.u. the ground state wave function is essentially composed of the P-func- 
_ tions Y', and ‘’,, whereas the wave function for the attractive, excited state is 
4 essentially composed of the S-functions ‘’, and V,. At R about 3 a.u. there is a 
_ strong mixing of different ‘Ys, and at larger R-values the ground state function 
is essentially composed of ‘; and ‘, and the attractive excited state of V; 
= and Y,. 


i ak 
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; 8. Results for the *y, states from superposition of configurations 


The calculations for the lowest 2X, state have been carried out in several steps, 
corresponding to the steps A, B and C used in Section 7 for the 7X; ground state. 
In Fig. 4 the result of step A, i.e. the calculation by use of ‘Y’; only, is given by the 
curve A of dashes. This curve is repulsive at large R-values, passes a maximum 
around R =3 a.u. and has a narrow minimum at Rk =1.5 a.u. The general character 
of the curve is preserved in the two following steps B and C, although the maximum 
becomes a little depressed and the minimum is somewhat deepened. 

The lowest 2X; state is the only excited state of Hz investigated by Dalgarno and 
McDowell [6]. Their most accurate treatment of this state is closely related to the 
present step A. They have only published their results of this treatment for R> 5 a.u. 
Within this range of R-values, containing only the repulsive part of the curve, the 
present results are in good agreement with those given in [6]. 

The result of step B, describing the interaction between the *X, states, derivable 
from the H- 48 and 3S states, where only 1s AO’s are involved, is given by the curve 
B of dots and dashes in Fig. 4. At large R-values this curve coincides with the the 
curve for step A, and at small R-values it coincides with the solid curve, representing 
the result of the following step C. This latter step, consisting of a superposition of 
W,, VY, and V,, describes the interaction between the 2; states, derivable from the 
H- 19, 38 and 8P states. In this approximation the minimum occurs at R = 1.54 a.u. 
and an energy value of —1.0546 a.u., corresponding to a D,-value of 0.0422 a.u. 
= 26.5 kcal/mole. é 

A step D, in which VP, VY, ‘V; and ¥’, are superposed, did not give any sensible 
result in the case of 2X¢ states. It is believed that this breakdown of the method 
depends on the fact that Y’, and Y’, separately give rise to energy values, which for 
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Fig. 5. Potential curves for the lowest ?X7, and *II,,., states of Hz together with the curve for 
the ground state of H,. 


R <3 a.u. only differ in the fifth decimal place. Now, it is estimated that the energy 
values are reliable only to the fourth decimal place. Moreover, in the superposition 
the rather unreliable difference between these two energy values is decisive. It is 
therefore not surprising that a step D does not work. 

As discussed in Section 7, several states of the same species can be calculated by 
use of the wave functions, given in (18). The final step of the calculation of 2X7 
states, i.e. the step C, gives thus rise to three different states. The energy values of 
all three states are listed in Table 5 of the Appendix. The two states, which are not 
represented in Fig. 4, are both repulsive. 


9. Results for the states *II,,¢, ‘Iu, and ‘Eu, 


The potential curves of the ?II states have been calculated partly by use of ¥',, of 
(20) and partly from a superposition of Y’, and ’,). The former case describes 21] 
states, derivable from the H- °P state. The latter case describes the interaction 
between II states, derivable from the H- !P and 8P states. 

The result for the lowest ?II, state is almost the same with and without superposi- 
tion of configurations. The discussion of the impossibility of superposing VY’; and Ys 
in a step D for the *X, case applies also to Y’, and ’,, in the 2IT, case. The energy 
values of all the II states are collected in Table 6 of the Appendix. One of the 2I1, 
states is repulsive. The other has a deep minimum at R = 1.64 a.u. with an energy 
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Fig. 6. Potential curves for the lowest *X7 , and ‘IT, , states of Hz together with the curve for 


the ground state of H,. 


value of —1.1050 a.u., only 0.0007 a. u. above the minimum of the 22; ground 
state. The corresponding D,-value is 0.1159 a.u. =72.7 kcal/mole. 

The lowest state in the best approximation of each of the types 7X, and 2II1,,, is 
shown in Fig. 5, together with the computed curve of the H, ground state. 

The *II,,,, potential curves are obtained by use of ’,, of (21). The computed energy 
values are listed in Table 7 of the Appendix. Both these states are found to be repul- 
sive. 

The *3,,, potential curves are computed in two different approximations: by use 
of ‘’, of (19) and from superposition of Y’, and ‘Y’,. The first case corresponds to the 
‘>; states, derivable from the 3S state of H~. The other case describes the interac- 
tion between the *D;., states, derivable from the 3S and °P states of H-. The pertinent 
energy values-are listed in Tables 8 and 9 of the Appendix. In both the 2, and the X, 
case one of the quartets is found to be repulsive. The other *X; state is found to have 
a very shallow minimum at R = 6.7 a.u. and an energy value of —0.9942 a.u., corre- 
sponding to a D,-value of 0.0051 a.u. = 3.2 kcal/mole. Because of this small value, 
comparable to van der Waals’ interaction, and because the whole potential curve 
of this state is found within the electron detachment continuum according to the 
present calculations, it cannot be claimed that the existence of this state has been 
ascertained. 

The lowest 4; state is found to have a minimum at Rk =3.4 a.u. and an energy 
value of —1.0164 a.u., corresponding to a D,-value of 0.0274 a.u. = 17.2 kcal/mole. 
This state seems therefore to correspond to an actual state of the Hz molecule. 
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- The lowest states of the species *D;,, and *II,,, are reproduced in Fig. 6 together 
with the computed curve of the H, ground state. 
10. Discussion of the results 


According to the results of the most accurate of the present calculations, presented 
in Figs. 5 and 6 and in Tables 3-9 of the Appendix several states of the molecule ion 


Hz should be stable towards dissociation into H and H-. Arranging these states. 


according to increasing value of their energy minima they are *X,, *II,, 2a aes 
De 

All the energy minima of the doublets are well below the energy value H,, of 
H(28) + H-(48), which is shown in Figs. 3-6 by the horizontal line, labelled H + H-. 
It is true that this line shows the value of #,,, computed according to Section 3, and 
is thus somewhat higher than the best known value of ,, [7], shown in Fig. 1 by 
the line H + Hy. It is however easily seen that all the energy minima of the doublets 
are also well below this best value of #,,. Therefore, no improvement of the wave 
functions could diminish the depths of these minima essentially, as did the improve- 
ments II, III, IV of the ground state function, summarized in Fig. 1. It seems there- 
fore to be reasonably certain that the states of Hz, listed above, actually correspond 
to descrete electronic energy levels of this molecule ion. However, the list may be 
incomplete and the given arrangement could be modified by future investigations, 
because the present calculations are particularly designed to give the best possible 
value of the ground state energy. 

It is interesting to compare this calculated energy level diagram of Hz with the 
corresponding diagram of H,. It is seen that the energy levels of Hz are very close 
together, approximately as the triplets of H,. This close spacing of the levels of Hz 
depends on the weak bonding of the extra electron, which in the orbital description is 
expressed by the diffusiveness of the outermost orbital. The close spacing of the 
energy levels is supposed to be a general feature of negative molecule ions. 

The potential curves for all discrete energy levels of Hg, calculated by the present 
author, refer to states obtained from interaction between H in its ground state and 
H_ in its ground state or in a low-lying excited state. As is seen from Fig. 5 all these 
curves lie comparatively close to the curve for the ground state of H,. Analogously, 
it may be suggested that states resulting from the interaction between H in an 
excited state and H~- in low-lying states will give rise to potential curves of Hz, 
lying close to the potential curve of H,, having one of the H atoms correspondingly 
excited. It is true that such hypothetical excited states of Hy should belong to the 
electron detachment continuum. Nevertheless, it is not definitely excluded that they 
might have a mean life, long enough to allow them to occur as intermediate steps in 
certain reactions, suggested below. 

For the discussion of the reactions, in which Hy could possibly participate, it is 
of interest to know (i) the dissociation energy D,, discussed above, (ii) the electron 
affinity of H,, EA (H,), (iii) the energy that must be supplied for vertical detachment 
of an electron from Hz, and (iv) the energy that must be supplied for vertical attach- 
ment of an electron to H,. The three latter quantities are determined by use of the 
potential curves for the ground states of both Hz and H,. Now, the potential curve 
of Hz is only known from the comparatively approximate theoretical calculations 
in [5] and [6] and in the present paper. As can be inferred from the corresponding 
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on Hy, performed so far, can by no means 


; estimate of the still unknown experimental potential. 


; with a potential curve of H,, calculated in a corresponding approximation, and 
not with the experimental potential curve. As discussed in Section 7 above, the 
potential curve of H, that corresponds to most of the treatments of Hz , carried out 
80 far, is the curve obtained from the simplest VB treatment, using 1s AO’s with 
: the constant value 1.00 of the orbital exponents. It is true that the orbital exponents 
have been varied both in [5] and [6] as well as in the present cases II, III and IV. 
_In view of the similarity between the results of the present cases III and A, the latter 
with the constant values (12) of the orbital exponents, it seems however to be justi- 
_ fied to compare all the different approximations of the potential curve of Hy with 
_ the above-mentioned H, curve. This is also the H, curve shown in Figs. 3-6. 
_ From a comparison of this potential curve of H, with the potential curve of Hz, 
_ arising from the present case ITI, the results of which are very close to the results 
_ of the best earlier investigation [6], the following values have been computed: 
_ EA(H,) = — 2.54 ev = —58.6 kcal/mole; vertical detachment energy of Hz = 
j + 0.46 ev = + 10.7 keai/mole; and vertical attachment energy of H, = + 8.87 ev = 
_ + 204.6 kcal/mole. It is seen that these quantities are widely different in the present 
_ case III. The same result can probably also be deduced from calculations according 
On the other hand, a comparison of the potential curve of H, with the best curve 
for the ground state of Hz, labelled D in Fig. 3, shows that in this case the three 
* above-mentioned quantities have very nearly the same absolute value, the attach- 
ment energy being positive and the other two negative, EA(H,) = —0.28 ev = 
— 6.4 keal/mole. 
As has been stressed above, the present theoretical potential curve D of Hg is 
_ still comparatively far from the correct experimental curve. It seems, however, as 
if the present result of a negative vertical detachment energy of Hz should be correct 
in view of the fact that He~ has a negative value of this energy, which by use of 
formula (1) can be estimated to be —0.5 ev. Thus, Hz should not be a stable, mo- 
lecular ion. 

By use of the results for Hz, discussed at some length above, a tentative explana- 
tion of the experimental results of Khvostenko and Dukelsky [3], [4] will now be 
suggested. According to the first of these papers, [3], the formation of H~ ions from 
collisions between electrons and hydrogen molecules has been measured as a func- 
tion of the energy of the electrons. Below 5 ev the yield of H~ ions was negligible. 
At 7.2 ev a maximum was found, which was attributed to H~ ions, formed by colli- 
sions between electrons and water molecules, a low H,O impurity of the H, gas 
being unavoidable, because apparently H,O was formed by interaction of H, gas 
molecules and O atoms of the glass wall. Another rather sharp peak was found at 
14.5 ev. Finally, after 18 ev it was found that the yield of H™ increased linearly with 
the electron energy. The authors attributed this last increase to a non-capture colli- 
sion, followed by dissociation of H, into H+ and H-. Moreover, they attributed the 
peak at 14.5 ev to resonance capture of electrons by hydrogen molecules, Hz being 
formed as an intermediate step. This suggestion seems to be confirmed by the later 
publication [4], according to which the authors obtained not only H- ions but also 
a low yield of Hz ions in the mass spectrometer, when the experiment was carried 


out under certain conditions. 


as 
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- These experiments are very difficult to understand in view of the earlier, theoreti- 
cal investigation of Hz [6] and the present case III, predicting a vertical electron 
attachment energy of H, as high as +8.9 ev and a dissociation energy of pity as 
low as 0.16 ev. The probability of a three-body collision, necessary to stabilize the 
Hz ion, seems to be too low to explain the measured amount of Hz. 

In view of the best approximation of the present investigation, predicting a 
vertical attachment energy of H, as small as + 0.28 ev and a dissociation energy 
of Hz as large as 2.5 ev, it does not seem to be too improbable that an electron 
could be captured on the collision with a hydrogen molecule, resulting in the forma- 
tion of an Hz ion. The present theoretical results are thus supported by the experi- 
ments reported in [3] and [4]. There is however one feature of the experiments, 
which is not explained by the assumption that the capture collision should give rise 
to the formation of Hz in its ground state, viz. the fact that the formation is only 
obtained within a rather narrow energy range around 14.5 ev. It is difficult to give 
an unambiguous explanation of this fact with the present, rather limited knowledge 
of the possible states of Hz. The following is only a tentative attempt. 

As is well known, electrons, having an energy within the range 11.8—15.6 ev, will 
excite H, molecules from the ground state to different, excited states, beginning 
with the triplet state a °X;. Electrons of an energy somewhat higher than 11.8 ev 
can therefore excite H, either to an excited, vibrational level of this triplet or to 
any vibrational level of a higher, electronic state. Now, according to the discussion 
above, it is not unlikely that there are excited energy levels of Hz close to the excited 
states of H,. The excited H, molecule could therefore possibly capture an electron 
in transition to an excited Hy state. If this is a 7I], or a 22> state, it can radiate in 
giving the 7X, ground state of Hz. Moreover, if the resulting electronic ground state 
were vibrationally excited, it might afterwards dissociate, giving H and H-. 

This tentative reaction scheme should also apply to the possible formation of 
H- by way of Hz in the solar corona [2]. These reactions will, however, be investi- 
gated further and discussed more closely in a forthcoming paper. 
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n: aod gives rise to integrals over 1s and 2p AO’s. The 
integrals were in this case computed on Besk, using a pro- 
written by Flodmark [12]. The exchange integrals of 1s AO’s only and with 
stant values (12) of the orbital exponents, as well as the hybrid and exchange 
rals, containing 2p AO’s, were obtained by courtesy of Dr. B. Ransil from com- 
utation on electronic digital computers, using the programmes written by members 
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Calculated energy values 


‘The energy > Fen of Hz, calculated for different internuclear distances and ac- 
- cording to different approximations, are collected in Tables 2-9. Values of inter- 
nuclear distances R are given in a.u., 1 a.u. = 0.52917 x 10-8 cm. Values of total 
energy, E are given in a.u., | a.u. = 97, 2097 ev = 627.54 kcal/mole. Dissociation 
energies D, = E,, — E are given both in a.u. and in kcal/mole. 

The values of R and # at the minimum were obtained by fitting a Morse curve to 
the three nearest points of the H(R) curve. The minimum values are therefore less 
accurate than the others, but are nevertheless given to the same decimal place in 
order to facilitate the comparison between the different approximations. 
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Table 2. Energy values FH (a.u.) for different internuclear distances R (a.u.) of the 
2); ground state of Hz according to the four cases, I, II, II, IV, discussed in 


Lov 


— Section 6. 

a “3 

‘ Energy values E of Hy according to case 

4 R 

Z I II III IV 

ye 

2 2.0 — 0.8528 — 0.9120 a — 
3.0 — 0.9108 — 0.9709 = = 

y 4.0 — 0.9134 —0.9719 —1.0141 Sis 

: 5.0 - — — —1.0191 

4 5.3 ae = = — 1.0192 
6.0 — 0.8899 — 0.9550 — 1.0185 — 1.0190 
8.0 — 0.8794 — 0.9464 — 1.0167 ST nile} 
10.0 — 0.8758 — 0.9421 — 1.0150 —1.0159 
oo — 0.8750 — 0.9404 —~ 1.0126 SLs 


ee 


Minimum: 


SN ae Me AS 


R (a.u.) "53156 3.47 5.72 5.42 

‘ E (a.u.) — 0.9152 — 0.9743 — 1.0185 — 1.0192 

az D, (a.u.) 0.0402 0.0339 0.0059 0.0059 
IBY (wank iaicie| 25.2 21.3 Be7 3.7 
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~1.0124 


1.64 1.65 
— 1.0518 — 1.1043 —l LOD 
2 (a.uU. 0.0394 0.0919 0.0933 
D, (kcal/mole) é 24.7 57.7 58.6 


Table 4. Energy values E (a.u.) for different internuclear distances R (a.u.) of the 
three excited ?X, states (1), (2), (3) of Hz, obtained from superposition of configu- 
rations according to the case D, discussed in Section 7. 


Energy values # of Ho for state 


R 

(1) (2) (3) 

1.0 — 0.8990 — 0.2518 + 2.5901 

1.5 — 1.0386 — 0.6813 — 0.4134 

2.0 = 1.0402 — 0.8439 — 0.6829 

3.0 — 1.0134 — 0.9314 — 0.8382 

4.0 — 0.9981 — 0.9413 — 0.8660 

6.0 — 0.9904 — 0.9539 — 0.8696 

8.0 — 0.9903 — 0.9621 — 0.8748 

fore) — 0.9891 — 0.9888 — 0.9217 

ce a 
Minimum: 

R (a.u.) 1.73 = — 
# (a.u.) = b0472 = — 
D, (a.u.) 0.0581 = — 


D, (kcal/mole) 36.5 = ae 
| DNS ERE 
; 


or 
bo 
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Table b. Energy values H (a.u.) for different internuclear distances R (a.u.) of the 
three X75 states (1), (2), (3), of Hz, obtained from superposition of configurations 
according to the cases A, B, C, discussed in Section 8. 


Energy values # of Hz , 
ee ree ee ee ee 


R For the lowest 4 state (1), According to case C for 
according to case the higher 4 states 

A B C (2) (3) 

1.0 — 0.9431 — 0.9549 — 0.9551 — 0.2942 — 0.1603 

1.5 — 1.0216 — 1.0537 — 1.0543 — 0.6727 — 0.5518 

2.0 — 0.9875 — 1.0316 — 1.0323 — 0.8369 — 0.7381 

3.0 — 0.9470 — 0.9606 — 0.9619 — 0.9519 — 0.8840 

4.0 — 0.9629 — 0.9637 — 0.9831 — 0.9583 — 0.9027 

6.0 — 0.9932 — 0.9943 — 0.9969 — 0.9845 — 0.8984 

8.0 —1.0041 — 1.0044 — 1.0048 — 0.9879 — 0.9069 
10.0 — 1.9081 — 1.0081 = = = 

co — 1.0124 — 1.0124 — 1.0124 — 0.9891 — 0.9217 

Minimum: 

R (a.u.) 1.5 1.54 1.54 at = 
E (a.u.) — 1.0216 — 1.0540 — 1.0546 — r= 
D, (a.u.) 0.0092 — 0.0416 0.0422 =e = 
D, (keal/mole) 5.8 26.1 26.5 = = 


Table 6. Energy values # (a.u.) for different internuclear distances R (a.u.) of the 
two II, states and the JI, state of Hz, discussed in Section 9. 


Energy values # of Hz for 


R ; *T1_, states from . 
Lowest *II,, superposition of Y, and ¥’,, | Lowest *I, 
state from state from 


Fro Lowest state | Higher state Pio 
1.0 — 0.9827 — 0.9849 + 3.4408 — 0.2913 
1.5 — 1.1006 — 1.1018 — 0.5828 — 0.6713 
2.0 — 1.0908 — 1.0927 — 0.8112 — 0.8358 
3.0 — 1.0287 — 1.0306 — 0.9347 — 0.9518 
4.0 — 0.9991 — 1.0005 — 0.9652 — 0.9803 
6.0 — 0.9887 — 0.9901 — 0.9757 — 0.9888 
8.0 — 0.9880 — 0.9898 — 0.9776 — 0.9893 
co — 0.9891 — 0.9891 — 0.9888 — 0.9891 
a eee ee an 
Minimum: | 

R (a.u. : 1.64 1.64 = = 

# (a.u.) —1,1036 — 1.1050 a = 

ve (a.u.) 0.1145 0.1159 = = 

. (kcal/mole) 71.9 7120 ee a ees es a mene = | 
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Table 8. Energy values E (a.u.) for different internuclear distances R (a.u.) of the 
two 42, states of Hz, discussed in Section 9. 


Energy values E of Hg for 


R States from super- 
Lowest state position of ‘Y", and ’, 
from Y’, 


Lowest state | Higher state 


1.0 — 0.2657 — 0.2845 + 1.9993 
2.0 —0.7905 — 0.9688 — 0.7694 
3.0 — 0.8897 — 1.0145 — 0.8830 
4.0 — 0.9085 — 1.0139 — 0.9067 
6.0 — 0.9176 — 0.9943 — 0.9176 
8.0 i — 0.9263 — 0.9857 — 0.9251 
co — 0.9217 — 0.9891 = O2N 
et ee ee Ee eee 
Minimum: 

R (a.u.) —= 3.4 = 

E (a.u.) = — 1.0164 = 

D, (a.u.) = 0.0274 = 

D, (keal/mole) a 17.2 A 


a 
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0.6902 
Se = 8248 
— 0.8592 
— 0.8637 
— 0.8635 | 
— 0.9217 


R (a.u.) 
E (a.u.) 
D, (a.a.) 
| D ama 
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